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Fault diagnoses of steam turbine using the exponential similarity 
measure of neutrosophic numbers 



Jun Ye* 

Department of Electrical and Information Engineering, Shaoxing University, 508 Huancheng West 
Road, Shaoxing, Zhejiang Province 312000, P.R. China 



Abstract 

Since the neutrosophic number consists of its determinate part d and its indeterminate part el 
denoted by N = d + el, it is very suitable for dealing with real problems with indeterminacy. 
Therefore, this paper proposes the exponential similarity measure of neutrosophic numbers and a 
fault diagnosis method of steam turbine by using the exponential similarity measure of neutrosophic 
numbers. By the similarity measure between the fault diagnosis patterns and a testing sample with 
neutrosophic numbers and its relation indices, we can determine the fault type and indicate fault 
trends. Then, the vibration fault diagnosis results of steam turbine demonstrate the effectiveness and 
rationality of the proposed diagnosis method. The proposed diagnosis method not only gives the 
main fault types of steam turbine but also provides useful information for multi-fault analyses and 
future fault trends. The developed diagnosis method is effective and reasonable in the fault diagnosis 
of steam turbine under an indeterminate environment. 

Keywords: Neutrosophic number; Exponential Similarity measure; Steam turbine; Fault diagnosis 



1. Introduction 

For the complex structure of steam turbine -generator sets, if there is a fault of the equipment, it 
will produce a chain reaction and cause the fault of other parts or equipments, and then will seriously 
impact the reliability of power generation. Therefore, one is convinced of the importance of fault 
diagnoses. The vibration fault of huge steam turbine -generator sets is usually a typical fault type. 
The fault diagnosis of the steam turbine realized by the frequency features is a simple and effective 
method. In many real situations, the testing data of the frequency features extracted from the 
vibration signals of steam turbine cannot provide deterministic values because the fault testing data 
obtained by experts are usually imprecise, uncertain and indeterminate due to a lack of data, time 
pressure, measurement errors, or the experts’ limited attention and knowledge. In some situations, 
the fault testing data usually contain the determinate information and the indeterminate information. 
Therefore, some researchers have developed the fault diagnosis methods of steam turbine based on 
the similarity measures of vague sets [1, 2, 3] and the cross entropy of vague sets [4] under an 
imprecise and uncertain environment. However, these diagnosis methods cannot deal with fault 
diagnosis problems with indeterminacy, which exist in real world, because the vague set cannot 
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include indeterminate and inconsistent information [5]. 

To deal with incomplete and indeterminate information in real world, Smarandache [6] firstly 
proposed the concept of a neutrosophic number, which is a subclass of neutrosophy. The 
neutrosophic number can be represented as A = d + el, which consists of its determinate part d and 
its indeterminate part el. In the worst scenario, A only contains its indeterminate part without its 
determinate part d, i.e. A = el. In the best scenario, there is only its determinate part d without its 
indeterminate part el, i.e. A = d. Obviously, it is very suitable for the expression of problems with 
determinate and indeterminate information. Therefore, the neutrosophic numbers can effectively 
represent the fault data with incomplete and indeterminate information. Recently, Ye [7] firstly 
provided a neutrosophic number tool for group decision-making problems with indeterminate 
information under a neutrosophic number environment and developed a de-neutrosophication 
method and a possibility degree ranking method for neutrosophic numbers from the probability 
viewpoint as a methodological support for group decision-making problems. Kang et al. [8] 
proposed the cosine similarity measure of neutrosophic numbers for misfire fault diagnosis of 
gasoline engines. Although the neutrosophic numbers have been defined in neutrosophic probability 
since 1996 [6], till now, little progress has been make for handling indeterminate problems by 
neutrosophic numbers in scientific and engineering areas. Furthermore, existing fault diagnosis 
methods for steam turbine cannot deal with fault problems with indeterminacy. In order to break 
through the predicament in engineering applications of neutrosophic numbers, this paper proposes an 
exponential similarity measure between neutrosophic numbers and a fault diagnosis method based 
on the exponential similarity measure for handling the fault diagnosis problems of steam turbine 
under a neutrosophic number environment. 

The remainder of this paper is organized as follows. Section 2 briefly describes some basic 
concepts of neutrosophic numbers. Section 3 presents the exponential similarity measure of 
neutrosophic numbers and investigates its properties. In Section 4, based on the exponential 
similarity measure of neutrosophic numbers, we establish a fault diagnosis method for the vibration 
fault diagnosis of steam turbine under an indeterminate environment and demonstrate the 
effectiveness and nationality of the proposed diagnosis method. Section 5 contains conclusions and 
future research. 

2. Neutrosophic numbers and their basic operational relations 

A neutrosophic number proposed by Smarandache [6] consists of the determinate part and the 
indeterminate part, which is denoted by N = d + el, where A and e are real numbers, and 7 is 
indeterminacy, such that P = I, 0x7 = 0 and V\ = undefined. 

For example, assume that there is a neutrosophic number A = 6 + 27. If 7 e [0, 0.4], it is 

equivalent to A e [6, 6.8] for sure A > 6, this means that the determinate part of A is 6, while the 

indeterminate part of A is 27 for 7 e [0, 0.4], which means the possibility for the number “A” to be 
within the interval [6, 6.8]. 

Let N\ — d\+ e\I and Ni = di + eil be two neutrosophic numbers for d\, e\, di, ei e R (R is all 
real numbers). Smarandache [6, 9, 10] introduced the operational relations of neutrosophic numbers 
as follows: 

(1) Ni + Ni = di + d 2 + (ei + ^2)7; 

(2) Ai — AO = di — di + (ei — ei)!',', 

(3) Ai x A 2 = did2 + {d\d2 + eidi + e\e2)P 
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(4) Nr = (d i + e\ I) 1 = dr + (2d\e\ + er)I\ 
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N, d,+e,I d, d^e, - d,e 2 T 

(5) — L = = — L + — ^ I for d 2 * 0 and d 2 *~ e 2 - 

N 2 d 2 +e 2 I d 2 d 2 (d 2 +e 2 ) 



( 6 ) yf~N\ - y[cl^+~ej - \ 



\yfd~i ~ (y[di + Jd i + e x )l 

d \ — d^ yj d i 1 ('. ) / 

— r/| + r/| + 4-^)7 

- - yjd^ e^I 



Let N\ — d\ + e\I and No = d 2 + e 2 I be two neutrosophic numbers. If d\, e\, d 2 , e 2 > 0, AO = d\ + 
e\l and N 2 = d 2 + e 2 I are called positive neutrosophic numbers. Then the relations between two 
positive neutrosophic numbers N\ and No are defined as follows: 

(1) N\ c No if and only if d\ + inf (e\I) < d 2 + inf (e 2 I) and d\ + sup (e\I) < d 2 + SU p (eol): 

(2) N i = No if and only if d\ + inf (e\I) = d 2 + inf (e 2 I) and d\ + sup (e\I) = d 2 + sup (e 2 I). 



3. Exponential similarity measure between neutrosophic numbers 

In this section, we propose the similarity measure between neutrosophic numbers based on 
exponential function. 

Based on an exponential function e~ x , we can give the definition of the exponential similarity 
measure between neutrosophic numbers. 

Definition 1. Let A = {Nai, Nai, ..., Na„ } and B = {Nb\, Nbi, ■■■, Nb,i} be two sets of neutrosophic 
numbers, where Naj = + ZajI and Nbj = dsj + esjl (j = 1,2, ...,«) for dAj, eAj, dBj, esj > 0. Then, the 

exponential similarity measure between A and B is defined as 



E(A,B) = -^e 



\\ d Aj +infi; e Aj I )-d Bj -inf( e Bj I ) + d Aj +sup (e Aj I )-d Bj -sup(^/ ) ] 



(i) 



Then, the exponential similarity measure satisfies the following proposition. 

Proposition 1. The defined exponential similarity measure E(A, B ) satisfies the following properties: 

(1) 0 < E(A, B) < 1; 

(2) E(A, B) = 1 if and only if A = B\ 

(3) E(A, B) = E(B, A); 

(4) If AcScC for a neutrosophic number C , then E(A, C) < E(A, B) and E(A, C) < E(B, C ). 

Proof: 

(1) It is obvious that the property is true. 

(2) If A = B, then there are Naj = dAj + eAjI = NBj = dBj + eBjl (j = 1,2, ..., n) for dAj, eAj, dBj, eBj > 0 and 
dAj = dBj and eAj = eBj, i.e. dAj + inf ( eAjI) = dBj + inf ( eBjl) and dAj + sup ( eAjI) = dBj + sup (eBjl). 
Hence there is E(A, B) = 1. 

If E(A. B) = 1 , there is dAj + inf (eAjI) = dBj + inf ( eBjl) and dAj + sup (eAjI) = dBj + sup (eBjl) since 
e° = 1 , i.e. dAj = dBj and eAj = eBj for / = 1,2, Obviously, there is A = B. 

(3) It is obvious that the property is true. 

(4) IfAcflc C, then this implies Eaj + inf (eAjI) ci dnj + inf (eBjl) and dAj + sup (eAjI) ci dBj + sup 
(eBjl). Hence, there are the following inequations: 
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I d Aj + inff e Aj I) - d Bj - inf( e Bj I)\ < \d Aj + inf( e Aj I) - d Cj - inf( e Cj I )\ , 

| d Aj + sup (e Aj I) - d Bj - sup f(e Bj I ) | < | d Aj + sup (e Aj I) - d Cj - sup(e Cj I)\, 

| d Bj + inf( e Bj I) - d Cj - inff e Cj I)\ < | d Aj + inff e Aj I) - d Cj - inf( e Cj I ) | , 

| d Bj + supO BJ I) - d Cj - sup(e Q /)| < | d Aj + sup (e AJ I) - d Cj - sup(e Q /)| . 

Thus, there are Ed A , C) < E(A, B) and E(A, C ) < E(B, C ) since the exponential function e~ x is a 
decreasing function for x e [0, oo). 

Therefore, the proofs of these properties are completed. |J§' : : 

In practical applications, the elements in the set of neutrosophic numbers are considered as 
different weights. Assume that the weight vector of the elements is W = (yvi, w’ 2 , ..., w„) T , w t e [0, 1] 

and vv = 1 • Thus, we can introduce the following weighted exponential similarity measure 
between A and B: 



E W (A,B) = X 

j = i 



w i e 



^ d A j+inKe Aj I)-d B j-m$ e Bj l ) + d^+supte^/ )-4 S j-sup(e s /)| 



(2) 



Obviously, the weighted similarity measure also satisfies the above properties (l)-(4) in 
Proposition 1. 



4. Vibration fault diagnosis of steam turbine using the exponential similarity measure 

In this section, a fault diagnosis method for the vibration fault diagnosis of steam turbine is 
proposed based on the exponential similarity measure of neutrosophic numbers. 

4.1 Fault diagnosis method based on the exponential similarity measure between neutrosophic 
numbers 

For a fault diagnosis problem with indeterminacy, neutrosophic numbers are most suitable for the 
expression of the indeterminate problem. Assume that a set of m fault patterns (fault knowledge) is P 
= {Pi, Pi, ..., Pm} and a set of n characteristics (attributes) is A = {Ai, A 2 , ..., A,,}. Then the 
information of a fault pattern Pk{k= 1,2, ..., m) with respect to a characteristic A, (/' = 1,2, ..., n ) is 
represented by a neutrosophic number Pk = {dn + e«/, dn + Ckil, ..., dk„ + ekd] ( k = 1, 2, ..., m), 
where dkj, etj > 0 (k = 1, 2, ..., m; j = 1, 2, ..., n). Then, a testing sample is represented by a 
neutrosophic number P, = [dn + e t \I, d t i + e r il, ..., d„, + e„J}, where d, h e,, > 0 (/ = 1,2, . .., n). 

The exponential measure value Vk{k= 1,2, ..., m) can be obtained by the following exponential 
similarity measure between P t and Pp. 



= E w (P k ,P t ) = £ 



Wj e 



\\d Kj +mi.e lj I)-d tJ -mXe ( jI )|+|</ t j+sup(ei,-7)-^-sup(c tf /)|J 



( 3 ) 



For convenient diagnosis, the measure values of Vk {k = 1,2, ..., m) are normalized into the 
values of relation indices within the interval [—1, 1] by 



2v, — v ■ — v 

k min max 



(4) 
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where v mx = max { v k } , v mn = min { v k } and r k e [- 1 , 1 ] for k - 1 , 2 , . . m. 

1 <k<m 1 <k<m 

Then, we can rank the relation indices in a decreasing order, determine the fault type, and 
indicate possible fault trends for the tested equipment. If the maximum value of the relation indices 
is rk = 1 , then we can determine that the testing sample P, should belong to the fault pattern Pi.. 

4.2. Vibration fault diagnosis of steam turbine 

Here, we introduce the vibration fault diagnosis of steam turbine by use of the proposed fault 
diagnosis method based on the exponential similarity measure of neutrosophic numbers to 
demonstrate the effectiveness of the proposed method. 

In the fault diagnosis problem of steam turbine [4], we consider a set of ten fault patterns P — 

{ f J i (Unbalance), /^(Pneumatic force couple), PAOffset center), TVOil-membrane oscillation), 
TVRadial impact friction of rotor), /^(Symbiosis looseness), /V Damage of antithrust bearing), 
TVS urge), TV Looseness of bearing block), Pn/Non-uniform bearing stiffness)} as the fault 
knowledge and a set of nine frequency ranges for different frequency spectrum A = [Ai(0.01-0.39/), 
A 2 (0.4-0.49/), A 3 (0.5/), A 4 (0.5 1-0.99/), A 5 (/), A 6 (2f), A 7 (3-5/), A 8 (Odd times of /), A 9 (High 
frequency > 5/)} as a characteristic set (attribute set) under operating frequency /. Then, the 
information of the fault knowledge Pu(k = 1,2, . . ., 10) with respect to the frequency range (attribute) 
A, (i = 1,2, ..., 9) can be introduced from [4], which are shown in Table 1. Then, assume that the 
weight of each characteristic A, is w, = 1/9 for i = 1,2, ..., 9. Each basic element (vague value [tkj, 

1 -fkj]) in a vague set Pa (k = 1,2, ..., 1 0; j — 1, 2, ..., 9) in Table 1 can be also considered as an 
interval value, then we can transform it into a neutrosophic number when indeterminacy I is within 
the interval [0, 0.02], as shown in Table 2. 

Table 1 Fault knowledge with vague values [4] 



Frequency range (f: operating frequency) 



Pk 

(Fault 

knowledge) 


Ai 

(0.01- 

0.39/) 


a 2 

(0.4- 

0.49/) 


a 3 

(0.5/) 


A 4 

(0.51- 

0.99/) 


a 5 

(/) 


Ae 

(2/) 


A 7 

(3-5/) 


As 

(Odd 

times 

of f) 


A 9 

(High 

frequency 

>5f) 


Pi 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.85, 


[0.04, 


[0.04, 


[0.00, 


[0.00, 


(Unbalance) 


0.00] 


0.00] 


0.00] 


0.00] 


1.00] 


0.06] 


0.07] 


0.00] 


0.00] 


P 2 (Pneumatic 


[0.00, 


[0.28, 


[0.09, 


[0.55, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.08, 


force couple) 


0.00] 


0.31] 


0.12] 


0.70] 


0.00] 


0.00] 


0.00] 


0.00] 


0.13] 


P 3 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.30, 


[0.40, 


[0.08, 


[0.00, 


[0.00, 


(Offset center) 


0.00] 


0.00] 


0.00] 


0.00] 


0.58] 


0.62] 


0.13] 


0.00] 


0.00] 


Pa 

(Oil-membran 


[0.09, 

A 1 1 1 


[0.78, 

A OAT 


[0.00, 
A A AT 


[0.08, 
a m 


[0.00, 
A A AT 


[0.00, 
A A AT 


[0.00, 
A A AT 


[0.00, 
A A AT 


[0.00, 
A A AT 



0 . 11 ] 0.82] 0 . 00 ] 0 . 11 ] 0 . 00 ] 0 . 00 ] 0 . 00 ] 0 . 00 ] 0 . 00 ] 



e oscillation) 
P 5 



(Radial 


[0.09, 


[0.09, 


[0.08, 


[0.09, 


[0.18, 


[0.08, 


[0.08, 


[0.08, 


[0.08, 


impact friction 
of rotor) 


0 . 12 ] 


0 . 11 ] 


0 . 12 ] 


0 . 12 ] 


0 . 21 ] 


0.13] 


0.13] 


0 . 12 ] 


0 . 12 ] 


Pr, 


[ 0 . 00 , 


[ 0 . 00 , 


[ 0 . 00 , 


[ 0 . 00 , 


[0.18, 


[ 0 . 12 , 


[0.37, 


[ 0 . 00 , 


[ 0 . 22 , 


(Symbiosis 


0 . 00 ] 


0 . 00 ] 


0 . 00 ] 


0 . 00 ] 


0 . 22 ] 


0.17] 


0.45] 


0 . 00 ] 


0.28] 
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looseness) 

P-, 




















r 7 

(Damage of 


[0.00, 


[0.00, 


[ 0 . 08 , 


[0.86, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


antithrust 

bearing) 


0.00] 


0.00] 


0.12] 


0 . 93 ] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


Pi 


[0.00, 


[ 0 . 27 , 


[ 0 . 08 , 


[ 0 . 54 , 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


(Surge) 


0.00] 


0 . 32 ] 


0.12] 


0 . 62 ] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


Pg 

(Looseness of 


[ 0 . 85 , 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[ 0 . 08 , 


[0.00, 


bearing block) 
Pio 


0 . 93 ] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


0.12] 


0.00] 


(Non-uniform 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[ 0 . 77 , 


[ 0 . 19 , 


[0.00, 


[0.00, 


bearing 

stiffness) 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


0 . 83 ] 


0 . 23 ] 


0.00] 


0.00] 



Table 2 Fault knowledge with neutrosophic numbers 








Frequency range (/: operating frequency) 






Pk 

(Fault 

knowledge) 


A\ 

(0.01- 

0 . 39 /) 


a 2 

( 0 . 4 - 

0 . 49 /) 


A3 

( 0 . 5 /) 


A 4 

( 0 . 51 - 

0 . 99 /) 


a 5 

(f) 


a 6 
(2 If) 


A 7 

( 3 - 5 /) 


Ai 

(Odd 

times 

of f) 


Ag 

(High 

frequency 

> 5 /) 


Pi 

(Unbalance) 


0 


0 


0 


0 


0 . 85 + 

7 . 5 / 


0 . 04 + 

2/ 


0 . 04 + 

1 . 5 / 


0 


0 


P 2 (Pneumatic 
force couple) 


0 


0 . 28 + 

1 . 5 / 


0 . 09 + 

1 . 5 / 


0 . 55 + 

7 . 5 / 


0 


0 


0 


0 


0 . 08 + 2 . 5 / 


Pi 

(Offset center) 


0 


0 


0 


0 


0 . 30 + 

14 / 


0 . 40 + 

11/ 


0 . 08 + 

2 . 5 / 


0 


0 


P 4 

(Oil-membran 
e oscillation) 


0 . 09 + 

/ 


0 . 78 + 

2/ 


0 


0 . 08 + 

1 . 5 / 


0 


0 


0 


0 


0 


P 5 




















(Radial 

impact friction 


0 . 09 + 

1 . 5 / 


0 . 09 + 

/ 


0 . 08 + 

2/ 


0 . 09 + 

1 . 5 / 


0 . 18 + 

1 . 5 / 


0 . 08 + 

2 . 5 / 


0 . 08 , 

2 . 5 / 


0 . 08 + 

21 


0 . 08 + 2 / 


of rotor) 

Pe 

(Symbiosis 

looseness) 


0 


0 


0 


0 


0 . 18 + 

21 


0.12+ 

2 . 5 / 


0 . 37 + 

4 / 


0 


0 . 22 + 3 / 


Pi 




















(Damage of 
antithrust 


0 


0 


0 . 08 + 

2/ 


0.86+ 

3 . 5 / 


0 


0 


0 


0 


0 


bearing) 

Pi 

(Surge) 


0 


0 . 27 + 

2 . 5 / 


0 . 08 + 

2/ 


0 . 54 + 

4 / 


0 


0 


0 


0 


0 


Pg 


0 . 85 + 


0 


0 


0 


0 


0 


0 


0 . 08 + 


0 
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(Non-uniform 0.77+ 0.19+ 

0 0 0 0 0 0 0 

bearing 3 1 21 
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In the vibration fault diagnosis of steam turbine, two real-testing samples are introduced from 
[4] as the form of neutrosophic numbers: P,\ = {0, 0, 0.1, 0.9,0, 0, 0, 0, 0} and Pa = {0.39, 0.07, 0, 
0.06, 0, 0.13,0, 0, 0.35}. 

For the fault diagnoses of the testing samples P t i and Pa, the exponential similarity measures 
and relation indices between Pu (k = 1,2, ..., 10) and P,, ( i = 1, 2) are calculated by Eqs. (3) and (4), 
as shown in Table 3. 



Table 3 Results of the relation indices and fault diagnoses 

Relation indices (n) Fault 



Pu diagnosis 





Pi 


Pi 


Pi 


Pa 


Ps 


Pe 


Pi 


Ps 


P9 


P 10 


result 


p, 1 


-0.5970 


0.1972 


-0.9014 


-0.5174 


-0.7707 


- 1.0000 


1.0000 


0.3037 


-0.5729 


-0.6780 


Pi 


Ptl 


-0.4736 


0.4815 


- 1.0000 


-0.4583 


0.7359 


0.7144 


-0.5086 


-0.6737 


1.0000 


-0.4397 


Pi 



For the first real-testing sample P t i, we can see from Table 3 that the fault type of the turbine is 
Pi due to maximum relation index ( 1 .0000), which indicates that the vibration fault of the turbine is 
firstly resulted from the damage of antithrust bearing. Then, the fault types of surge (Ps) and 
pneumatic force couple (Pi) have possibility because the relation indices of the fault types Ps and Pi 
are positive. While the fault types of Pi, P 3 , P 4 , P 5 , P 6 , P<>, and P 10 have very low possibility due to 
the negative relation indices. By actual checking, we discover that one of antithrust bearings is 
damage. Therefore, it causes the violent vibration of the turbine. Flereby, the ranking order of all 
faults is P 7 — yP s — y P 2 — yP 4 — yP 9 — yP 1 — yP 10 — yP 5 — yP 3 — yP 6. 

For the second real-testing sample Pa, we can see from Tale 3 that the vibration fault of the 
turbine is firstly resulted from the looseness of bearing block (Pg), next the radial impact friction of 
rotor (P 5 ) and the symbiosis looseness (Pt,) have high possibility because their relation indices are 
more than 0.5, then pneumatic force couple (Pi) has possibility because the relation index of the 
fault type Pi is positive. Clearly, the fault types of Pi, P 3 , P 4 , Pi, Ps and Pin have very low possibility 
due to the negative relation indices. By actual checking, we discover the friction between the rotor 
and cylinder body in the turbine, and then the vibration values of four ground bolts of the bearing 
between the turbine and the gearbox are very different. We also discover that the gap between the 
nuts and the bearing block is oversize. Thus, the looseness of the bearing block causes the violent 
vibration of the turbine, which is in accordance with the actual fault. Flereby, the ranking order of all 
faults is P 9 — yP 5 — yP 6 — yP 2 — yP 1 0 — yP 4 — yP 1 — yP 1 — yP 8 — yP 3 - 

Obviously, the results of all fault diagnoses are the same as the actual faults of steam turbine. 
The fault diagnosis results of the turbine show that the proposed method not only indicates the main 
fault types of the turbine, but also provides useful information for multi-fault analyses and future 
fault trends. The proposed methods in this paper are effective and reasonable in the fault diagnoses. 
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5. Conclusion 

This paper proposed the exponential similarity measure between neutrosophic numbers and its 
fault diagnosis method for the vibration fault diagnosis of steam turbine under a neutrosophic 
number environment. The fault diagnosis results demonstrated the effectiveness and rationality of 
the proposed method. The proposed diagnosis method can not only determinate the main fault type 
of steam turbine but also predict future fault trends according to the relation indices. The proposed 
diagnosis method can handle fault diagnosis problems with indeterminacy and is simpler and easier 
than existing diagnosis methods, while the existing diagnosis methods cannot deal with fault 
diagnosis problems with indeterminacy. The diagnosis method proposed in this paper extends 
existing diagnosis methods and provides a new way for fault diagnoses under an indeterminate 
environment. In the future, the developed diagnosis method will be extended to other fault diagnoses, 
such as vibration faults of aircraft engines and gearboxes. 
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Table(l) 

Click here to download Table: nnesm_stfdl_tablel.doc 



Table 1 Fault knowledge with vague values [ 4 ] 



Frequency range (/: operating frequency) 



Pk 

(Fault 

knowledge) 


A\ 

(0.01- 

0 . 39 /) 


a 2 

( 0 . 4 - 

0 . 49 /) 


A3 

( 0 . 5 /) 


A4 

( 0 . 51 - 

0 . 99 /) 


a 5 

(/) 


A 6 
(2 f) 


A 7 

( 3 - 5 /) 


As 

(Odd 

times 

of f) 


Ag 

(High 

frequency 

> 5 /) 


Pi 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[ 0 . 85 , 


[ 0 . 04 , 


[ 0 . 04 , 


[0.00, 


[0.00, 


(Unbalance) 


0.00] 


0.00] 


0.00] 


0.00] 


1.00] 


0 . 06 ] 


0 . 07 ] 


0.00] 


0.00] 


F2 (Pneumatic 


[0.00, 


[ 0 . 28 , 


[ 0 . 09 , 


[ 0 . 55 , 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[ 0 . 08 , 


force couple) 


0.00] 


0 . 31 ] 


0.12] 


0 . 70 ] 


0.00] 


0.00] 


0.00] 


0.00] 


0 . 13 ] 


Pi 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[ 0 . 30 , 


[ 0 . 40 , 


[ 0 . 08 , 


[0.00, 


[0.00, 


(Offset center) 


0.00] 


0.00] 


0.00] 


0.00] 


0 . 58 ] 


0 . 62 ] 


0 . 13 ] 


0.00] 


0.00] 


P 4 

(Oil-membran 


[ 0 . 09 , 


[ 0 . 78 , 


[0.00, 


[ 0 . 08 , 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


e oscillation) 

P5 


0.11] 


0 . 82 ] 


0.00] 


0.11] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


(Radial 


[ 0 . 09 , 


[ 0 . 09 , 


[ 0 . 08 , 


[ 0 . 09 , 


[ 0 . 18 , 


[ 0 . 08 , 


[ 0 . 08 , 


[ 0 . 08 , 


[ 0 . 08 , 


impact friction 
of rotor) 


0.12] 


0.11] 


0.12] 


0.12] 


0.21] 


0 . 13 ] 


0 . 13 ] 


0.12] 


0.12] 


Pb 

(Symbiosis 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[ 0 . 18 , 


[0.12, 


[ 0 . 37 , 


[0.00, 


[0.22, 


looseness) 

Pi 


0.00] 


0.00] 


0.00] 


0.00] 


0.22] 


0 . 17 ] 


0 . 45 ] 


0.00] 


0 . 28 ] 


n 

(Damage of 


[0.00, 


[0.00, 


[ 0 . 08 , 


[0.86, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


antithrust 

bearing) 


0.00] 


0.00] 


0.12] 


0 . 93 ] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


Pi 


[0.00, 


[ 0 . 27 , 


[ 0 . 08 , 


[ 0 . 54 , 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


(Surge) 


0.00] 


0 . 32 ] 


0.12] 


0 . 62 ] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


P9 

(Looseness of 


[ 0 . 85 , 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[ 0 . 08 , 


[0.00, 


bearing block) 

Pio 


0 . 93 ] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


0.12] 


0.00] 


(Non-uniform 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[0.00, 


[ 0 . 77 , 


[ 0 . 19 , 


[0.00, 


[0.00, 


bearing 


0.00] 


0.00] 


0.00] 


0.00] 


0.00] 


0 . 83 ] 


0 . 23 ] 


0.00] 


0.00] 



stiffness) 



Table(2) 

Click here to download Table: nnesm_stfd1_table2.doc 



Table 2 Fault knowledge with neutrosophic numbers 









Frequency range (f: operating frequency) 






Pk 

(Fault 

knowledge) 


A\ 

( 0 . 01 - 

0.39/) 


a 2 

(0.4- 

0.49/) 


A 3 

(0.5/) 


A 4 

(0.51- 

0.99/) 


a 5 

(/) 


A 6 

( 2 /) 


A 7 

(3-5/) 


As 

(Odd 

times 

of f) 


A 9 

(High 

frequency 

>5/) 


Pi 

(Unbalance) 


0 


0 


0 


0 


0.85+ 

7.5/ 


0.04+ 

2 / 


0.04+ 

1.5/ 


0 


0 


P2 (Pneumatic 
force couple) 


0 


0.28+ 

1.5/ 


0.09+ 

1.5/ 


0.55+ 

7.5/ 


0 


0 


0 


0 


0.08+2.5/ 


Pi 

(Offset center) 


0 


0 


0 


0 


0.30+ 

14/ 


0.40+ 

11 / 


0.08+ 

2.5/ 


0 


0 


Pa 

(Oil-membran 
e oscillation) 


0.09+ 

/ 


0.78+ 

2 / 


0 


0.08+ 

1.5/ 


0 


0 


0 


0 


0 


Ps 




















(Radial 

impact friction 


0.09+ 

1.5/ 


0.09+ 

/ 


0.08+ 

2 / 


0.09+ 

1.5/ 


0.18+ 

1.5/ 


0.08+ 

2.5/ 


0.08, 

2.5/ 


0.08+ 

21 


0.08+2/ 


of rotor) 
Pe 

(Symbiosis 

looseness) 


0 


0 


0 


0 


0.18+ 

2 / 


0 . 12 + 

2.5/ 


0.37+ 

4/ 


0 


0.22+3/ 


Pi 




















(Damage of 
antithrust 


0 


0 


0.08+ 

2 / 


0 . 86 + 

3.5/ 


0 


0 


0 


0 


0 


bearing) 

(Surge) 


0 


0.27+ 

2.5/ 


0.08+ 

2 / 


0.54+ 

4/ 


0 


0 


0 


0 


0 


P9 

(Looseness of 
bearing block) 


0.85+ 

4/ 


0 


0 


0 


0 


0 


0 


0.08+ 

21 


0 


Pio 




















(Non-uniform 

bearing 


0 


0 


0 


0 


0 


0.77+ 

3/ 


0.19+ 

2 / 


0 


0 


stiffness) 





















Table(3) 

Click here to download Table: nnesm_stfd1_table3.doc 



Table 3 Results of the relation indices and fault diagnoses 

Relation indices (rk) Fault 



Pu diagnosis 





Pi 


Pi 


Pi 


Pa 


Pi 


P6 


Pi 


Ps 


P 9 


P 10 


result 


Pa 


-0.5970 


0.1972 


-0.9014 


-0.5174 


-0.7707 


- 1.0000 


1.0000 


0.3037 


-0.5729 


-0.6780 


Pi 


Pa 


-0.4736 


0.4815 


- 1.0000 


-0.4583 


0.7359 


0.7144 


-0.5086 


-0.6737 


1.0000 


-0.4397 


P 9 



